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SCALING LIMIT RESULTS FOR THE SUM OF MANY INVERSE 

LEVY SUBORDINATORS 

INGEMAR KAJ AND ANDERS MARTIN-LOF 



Abstract. The first passage time process of a Levy subordinator with heavy-tailed 
Levy measure has long-range dependent paths. The random fluctuations that appear 



(N 

o 

?H ' under two natural schemes of summation and time scaling of such stochastic processes 

are shown to converge weakly. The limit process is fractional Brownian motion in one 
case and a non-Gaussian and non-stable process in the other case. The latter appears 
to be of independent interest as a random process that arises under the influence 
C 3 ' of coexisting Gaussian and stable domains of attraction and is known from other 

£f) applications to provide a bridge between fractional Brownian motion and stable Levy 

motion. 

+-> 



1. Introduction and statement of results 



A Levy subordinator {X t ,t > 0} is a real- valued random process with independent 
and stationary increments and increasing pure-jump trajectories. The inverse process 
{T x , x > 0} defined by the first passage times T x = inf{£ > : X t > x} has nondecreasing, 
trajectories, where the lengths of the flat pieces of {T x } correspond to the jump sizes of 
{X t }. The dependence structure in the paths of the inverse process is entirely different 
from that of the Levy subordinator, since big jumps in the Levy process may cause strong 
dependencies that last over a considerable period of evolution of the path of its inverse. In 
\Q , this paper we take a scaling approach to study the nature of the random fluctuations that 

build up as a result of such long-memory effects. By superposing a large number of paths 
of the inverse Levy process and simultaneously scale the time parameter of the process, 
we obtain scaling limit results for the centered and normalized superposition process. 

In somewhat more detail, our starting point is a Levy subordinator with Levy measure 
v(dx) of regularly varying tail with index 1 + /3, < /3 < 1. In particular, /i := j xv{dx) < 
oo. The initial distribution of the subordinator process is chosen such that the resulting 
inverse process has stationary increments and expected value E{T X ) = x/fi. Letting 
{^x}«>i be a collection of independent copies of {T x } 7 our main result is the derivation of 
a limit process for the summation scheme 

m 1 

as both m and a — a m tend to infinity in such a way that m is of the same order of 
magnitude as a* 3 , modulo slowly varying functions. The reason for this choice of scaling is 
to attempt to trace the superposition process on a time scale that captures the size of the 
fluctuations around its mean. In the asymptotic limit appears a non-Gaussian, non-stable 
process with long-range dependence, which is known to arise also in other related models 
and has been called fractional Poisson motion, |5], [5], [13], [5]. The general study [5] of 
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higher-order moment measures for heavy-tailed renewal point processes, provides a unified 
framework of [5] and the present work. 

To give a heuristic context for the topics of interest in this work, let us recall the 
following limit result for Levy processes. Writing a = 1 + /3, the centered and scaled 
process (X t — (it) ft 1 / 01 converges in distribution as t — > oo to a random variable Z a , 
having a stable distribution with stable index a. If we write T x for the overshoot at x, so 
that Xt x = x + T X , then 

T x -x/n X Tx -iiT x (T x \V<*l T x 



m 



x 1 / 01 T x /a ^ x ) n fix 1 / 01 ' 

In this relation, T x /x — > 1/fi as x — > oo by the law of large numbers. It can be shown 
moreover that the second term on the right hand side is a remainder term with T x jx x l a — > 
as x — > oo. Therefore (T x — x/^/x 1 ^ 01 converges in distribution to —Z a /fi 1+1 / a as 
x — > oo. Proceeding heuristically, with m ~ or we may rewrite the superposition process 
cither as 

1 _ i i r ax 1 m 1 

- y2( T ax - -«#) — f-g^ / —fH y2( dT u - -du) 

a *—' a a L p l z In m L l z ■f— ' u 

8=1 ^ * /U 1 = 1 ^ 

or 

a Z^^ J - „ a ^ m l/(l+/9)Z^ a l/(l+/3) • 

Z=l 2=1 

The hrst representation emphasizes a sequence of random variables in the domain of 
attraction of a Gaussian law (m — > oo with a fixed) . The second representation highlights 
a sequence in the domain of attraction of a stable law with index 1 + /3 (a — >• oo with 
m fixed), which is the type of convergence just discussed above. For the limit regime of 
interest in our case Gaussian and stable attraction appear to coexist and both influence 
the resulting limit process. 

The main result (Theorem 2 below) is a scaling limit theorem for the intermediate 
type rescaling regime indicated above. In parallel to this we discuss the scaling regime 
of Gaussian predominance, leading to fractional Brownian motion in the limit (Theorem 
1). Scaling limit results with fractional Brownian fluctuations are known for a variety of 
models, such as modeling random variation in aggregated data traffic streams. For an 
introduction and overview of these topics and discussion of the modeling context, as well 
as detailed statements and derivations of such results, see e.g. [TBI fTTl [T3] . 

The model is introduced in detail and all results are stated in Section 1 of the paper. 
We then focus on the proof of Theorem 2 for the intermediate scaling regime, starting with 
the analysis of marginal distributions in Section 2. The main technique we use for the 
study of the one-dimensional distributions of the scaled processes and their limit behavior 
is that of double transforms in the sense of taking Laplace transforms in the time variable 
of the logarithmic moment generating function of the random variables. In Section 3 we 
continue with a study of the finite-dimensional distributions, which are obtained from 
recursive sets of integral equations for the finite-dimensional cumulant functions. Finally 
in Section 4 we provide a summary of the proof of Theorem 1 for Gaussian scaling, where 
each parallel step turns out to be simpler, and the proof of tightness. 

1.1. A Levy subordinator and its inverse. Let {X t ,t > 0}, Xq — 0, denote a Levy 
subordinator with right-continuous paths, having drift zero and Levy measure f(a, b) = 
J v(dx) with no atom at zero, such that 

/•oo />oo 

(1) / (1 A x)v{dx) < oo and \i = / xv(dx) < oo, 

Jo Jo 
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which implies that the first moment is finite, E(X t ) = fit < oo. The Laplace transform is 
given by — \nE(e~ uXt ) — t$(u), u > 0, with Laplace exponent 



§(tt)= / {l-e- ux )v{dx). 
Jo 



for which E(e~ uX °) = 




(3) 


£( e - uXi 



/() 

Let Xt — X n + X t denote the corresponding delayed subordinator process with general 
initial distribution Xq assumed to be independent of {X t }. We will study the case when 
Xq > has distribution function 

(2) P(X <x) = - / u(da)dy, 

and so 

f ) = — $(u) exp{-i$(u)l u > 0. 

flU 

Next we introduce the first passage process of the subordinator. Useful references are 
Bertoin pQ , [5] . Van Harn and Steutel, [TTJ , investigate stationarity properties of delayed 
subordinators and derive closely related results to those in Lemma 1 and Lemma 3 below. 
The entrance time of the Levy process {X t } into a set B is defined by Tb = inf{£ > : 
X t € B}. For any open set B, Tb is a stopping time. The first passage time T x = T( Xi00 ) 
strictly above a level x is the entrance time into (x, oo), that is 

Ta = inf{£ > : X t > x}, x>0, 

which is a right-continuous function with left limits. Since Xt f oo as t f oo, we have 
T x < oo for all x and P^ < i) = P(X t > a;), X, t > 0. Also, 

S(T X ) = / P{T X >t)dt = E l{x t < x ) dt. 

Jo Jo 



Hence 

/■oo />oo 

/ ue~ ux E{T x )dx= I E(e- uXt )dt = 

jJM 



/ ue- ux E{T x )dx= E(e- uXt )dt= — 

Jo Jo 



in view of ([3]), and therefore 

£(T X ) = -x. 

We call {T^} the inverse Levy subordinator and the process 

f x = inf{t > : X t > x}, x > 0, 

the pure inverse Levy subordinator. The path-regularity of {T x } is determined by the 
distribution of the small jumps of {Xt}, manifest in the asymptotic behavior of v{dx) 
for x close to 0. The nature of the paths varies considerably with compound Poisson 
processes as one extreme case. These are the subordinators for which the Levy measure 
is finite on the positive half line and the passage time process has piece-wise constant 
trajectories of lengths drawn from the probability distribution i/(0,x]/i/(0, oo), x > 0, and 
with exponentially distributed jumps. On the other hand, if the number a = sup{a > : 
lim^-^oo A -a< 5(A) = oo}, known as the lower index of the subordinator, is positive (and 
< 1 because of assumption ((JJ) then the inverse process {T x } is a.s. 7-Holder continuous 
on any compact interval, for each index 7 < a, see [JJ, Ch. 3 (Xq = 0). 

The scaling problem studied in this work involves weak convergence in the sense of 
convergence of finite-dimensional distributions plus a tightness property in the space C = 
C[0, 00) of continuous random processes. The path space is such that each C[0, T], T > 0, 
is equipped with the topology of convergence in supremum norm. 

We will prove below the following 



4 INGEMAR KAJ AND ANDERS MARTIN-LOF 

Lemma 1. The inverse subordinator process {T x ,x > 0} has stationary increments. 

1.2. Scaling limit theorem. Our basic assumption is that the Levy measure v is regu- 
larly varying at infinity with index 1 + /3, < /3 < 1, i.e. 

,oo 1 

(4) / v(dy) j—gi(x), a: -» oo, 

J a; ^ 

where L is a slowly varying function and we write f(x) ~ fl 1 ^) if / and g are positive 
functions and f(x)/g(x) — >• 1 as x — > oo. The summation schemes to be applied involve 
speeding up the time parameter using a rescaling sequence a m — ¥ oo, either such that 

(5) mL(a m )/a^ n — > oo, to — )■ oo, 
or such that 

(6) mL(a m ) / a„ — > <r /i, m — >• oo, 

where c, < c < oo, is an additional parameter that signifies the relative change of scales 
of size and time. In addition, we assume that the lower index of the Levy measure v 
satisfies 

(7) a = sup{a > : lim A~ Q $(A) = oo} > f3. 

A— >-oo 

Theorem 1. Under assumptions O) and Hj), let a m be a sequence such that a m — » oo as 
m — > oo and ^ holds, and define b m by 

(8) b 2 m = ma 2 ~ fi ! L(a m ) / '/z. 

Then, in the sense of weak convergence of random processes in C, 

( 9 ) \jT YfiL* ~ ~ amX ^ x ^°}^ {^^pB H {x), x > 0}, 
where 

and Bh is standard fractional Brownian motion with Hurst index H , i.e. the Gaussian pro- 
cess with stationary increments, variance V(Bjj(i)) = t 2H and continuous sample paths. 

n . n n 

iog£ex P { e e,B H ( Xl )} = j EE wtfr* + x T p (*< - ^) 2-/s ). 

i—1 i—1 j—1 

where = x$ < x\ < ■ ■ ■ < x n , n > 1. 

Theorem 2. Under the assumptions (01) and Fw, if a m is a sequence such that a m — > oo 
and (0|) holds for some constant c > as to — > oo 7 then 

m 

(10) {— E( T L. - ~ a ™x), x > °} =► {-M _1 c^(o;/c), x > 0}, 

in i/ie sense of weak convergence in C. Here {Yp(x),x > 0} is a zero mean stochastic 
process with continuous paths and finite- dimensional distributions characterized by the 
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cumulant generating function 



r 1 1 f f 



e^u-P dudv 



1 n — 1 n j—1 

^ Z=l J=t+1 & = 2+l 

(11) x/ / e 6 > u e eiV {xj-i -x. + u + vy 13 dudv, 



JO JO 

where = Xq < x\ < ■ ■ ■ < x n , and Axi — Xi — ajj-i, i = 1, . . . , n. 

Remarks, a) It is known that {Yp } has a simple representation as the stochastic integral 



Yp(x)= / / l{ s <; / < s+ ^}dyiV(ds,di)) 

Jj?,x[0,oc) JO 

where N(ds,dv) = N(ds,dv) — n(ds,dv) and N(ds,dv) is a Poisson random measure on 
R x [0, oo ) with intensity measure n(ds,dv) = (1 + (3)dsv~ 2 ~ l3 dv, [HI HI [13]. We have 
not been able to find however a method of proof of the present results which utilizes more 
directly this inherent Poisson structure of the model. 

b) The process {Y@ (x)} has been derived in [9] as a limit process in the setting of a su- 
perposition of independent renewal processes with stationary increments and heavy-tailed 
inter-renewal distribution, and in [12] and [13] for an infinite source Poisson process with 
heavy-tailed activity periods. The motivation is partly from modeling the total traffic load 
generated by many independent sources at an arrival point in a data traffic network. In 
these references condition ([5|) is called fast connection rate and ([6]) intermediate connection 
rate. They are compared to an alternative third scaling regime of slow connection rate, for 
which the limit process turns out to be a stable Levy process with stable index a = 1 + j3, 
see also et al. [14] or et al. [17]. 

c) Proofs of the following properties among others can be found in Gaigalas and Kaj 
[9]. The process {Yg} has stationary increments and continuous trajectories. The process 
is not self-similar. The higher moments are of the order E{Yp{x)) ~ const x k ~@ , k > 2, 
for large x. Specifically, the second-order properties (mean, variance, covariance) are the 
same (modulo constants) as those for fractional Brownian motion, whereas higher order 
moments are different. For example, {Yg} is positively skewed. The paths are 7-Holder 
continuous for all 7 < 1 — [3/2 (not 7 < 1 as claimed in [5]). 

d) The renewal processes studied in [9] can be viewed as discrete local time processes 
of discrete regenerative sets (ranges of compound Poisson subordinators). In this light, 
the present situation is the natural analogue for continuous local time processes of perfect 
regenerative sets (ranges of subordinators that are not compound Poisson) . One can expect 
the scaling limits to transfer since they are large-time asymptotics which should not depend 
on the local structure. Some relevant references for the connections of regenerative sets 
and subordinators are [7], and [ID] . 

2. Analysis of the marginal distrbution 

As a preliminary for the proof of Theorem 2 we observe the following properties of the 
functions introduced in (|llj) . which are straightforward to verify. 
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Lemma 2. Relation ill]) defines a consistent family of finite- dimensional distributions, 
such that for any c > 



n 

\ogE exp {^e^cYpixt/c) ~ cY fj (x^i/c))j 

n 



The main part of the proofs of Theorem 1 and Theorem 2 consists in establishing 
convergence of the scaled n-point cumulant functions 



log£exp { J2 Si^~ X>i^ " -VmXi)} 

t=l ° m fc=l V J 

n ft 1 

(12) = mE \ exp { ^ -± {T amXi - -a m Xi) } - l] + 0(l/m) 

toward the corresponding functionals of the limit processes. As a preparation we study 
the joint distribution (T X ,T X ), where {T^} is the overshoot process, and other properties 
of the one-dimensional marginal distributions of T x . 

2.1. Marginal distributions and the overshoot process. The overshoot process {T x , x > 
0} associated with the first passage time {T x } is defined for x > by 

T x = X Tm - x, 

and represents at time x the remaining time until the next point of increase of the inverse 
subordinator. The following Lemma is a special case of a result valid for general Levy 
processes adaptated to the case of a general initial distribution Xq. For a proof see 
Theorem 49.2 in Sato [IS]. 

Lemma 3. For u > 0, 9 < $(w), and v > with v ^ u, 

J u — v \ /j,v $(u) — 8 fiu J 

Proof of Lemma [I] In Lemma [3] take u > and v > 0, u ^ v, and let 9 = 0. We obtain 



f 

Jo 



ue-™E(e-^)dx- U W v) * (w) $( ^ - $(u) 



u — v \ fiv Q(u) [iu ) fj,v 



Hence, for any x > 0, T x = Xq, Consequently, for each x the increment process T x+y — T x , 
y > 0, begins with a flat period for a duration of time having the distribution Xq, which 
is just the same behavior as the original process T x , x > 0. To formalize the argument, 
note 

P(T x+y ~T x >t) = P(T X < y, T x+y -T x >t) = P(T X < y, X Ti+ t - X n < y - T x ). 

Since F^ = Xt x — x is independent of XT x +t — Xt x and Xr x +t — Xt x = -X* it follows that 
P(T x+y -T x >t) = P{X < y, X Tsc+t -X Ti <y- X ) = P(X t <y) = P(T y > t). 

□ 
Lemma 4. For u > and 8 < $(u), 

(13) / ue- ux E{e al *)dx = l 



JT^^_, & *(«) 



$(u) — 9 fiu 
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Also, for u > and 9 > —\iu, 

/>oo 

(14) / ue~ ux E{e e{T *~ x/tl) - 1) dx 

Jo 



(jiu + 9f 



/;?/ 



and 

(15) 



ue~ ux E(e 9{T *- x/ ^ -l)dx = - 

fiu + 9 



$(« + 9 /n) - 



$(u + 9/ll)-9 



1 



Proof. Relation (fj~3j) follows by letting «^0in Lemma [3] and using that <&(v)/v — > fi 
in this limit. 

The remaining calculations, involving the random variables T x — x/fi and T x — x/fi, 
follow from (p~3|) and the analogous expression 



UX EV„ST X 



ue~ ux E[e 



dx 



$(u) 



n 



, 9(u)-0' 

where we note cf>(u + 6/fi) < fiu + 9 for all 9 such that u + 9/fi > 0. 
Lemma 5. The function 

£( e *(T.-*//i) _ -^ ^ > Qj 

is nonnegative for any real parameter 9 and differentiable and nondecreasing with respect 
to the variable x. The derivative with respect to x is given by 

± E , e 6(T*-x/n) _ 1) = 6e -e*/» E ( e 6T x _ e er.y > 
dx 

Proof. The nonnegativity follows from Jensen's inequality. It follows from (|14l) . (fT5|) 
and the uniqueness property of Laplace transforms that E{e 6<yT:c ~ x I ^ — 1) is obtained as 
the convolution of E(e 9 ^ Tx ~ x ^ IJ ^ — 1) with the exponential e~ 9x ^. Hence 



E(( 



9(T m -x/ij) 



1 = - I , 

M Jo 



-^-v)/^E(e 9 ^- y /^-l)dy. 



The left hand side is differentiable in x with derivative 

± E t e e{T m -x/ti _ l) = _fL E ( e KT.-*M _ !) + i E f e e(.f m -x/ ll ) _ q 
dx jj, fx 



-e- ex /»E(e 6T * 



e eT *). 



Now we observe that the processes T x and T x can be constructed on the same probability 
space by a shift of size X so that T is a copy of T with the first point of increase in X 



rather than in 0. In particular P{T X > T x ) = 1. Hence 9E(e 



QT X d st x 



> for any 9. 



□ 



Lemma 6. For x > 0. 

i) x/n < E{f x ) < 



$(!/»)■ 



V(X, CO J 

iii) — Var(T x ) = -£(?„ - .T//i) > 0, 



dx 



2e 



H 



iv) Vax(T x ) < — / ^(1/y)- 1 dy - (x/^i) 2 . 
A* Jo 
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Proof. For i), it was noticed in the proof of Lemma [3] that the processes T x and f x 
could be constructed such that T x > T x almost surely. Hence E(T X ) > E(T X ) = xj [i. 
Moreover, 

/■OO /"OO _ /"OO 

E{T X ) = P{X t <x)dt< eE(e- Xt/x )dt = e e ~**( 1 /-) dt = e /$(l/x). 

Jo Jo Jo 

Inequality ii) follows from 

$(!/«) = - / e - u/x v(u, oo) du>- / e- u,x v{u, oo) du > !/(*, oo)(l - e" 1 ). 



To prove iii) and iv), differentiate twice with respect to 9 in (fT4| to obtain 

v ' Jo (/z U ) 2 v$( U ) y 

Similarly, using (fT5|) . 

/°° W e-«^(f x - x/a«) c*r = — (-j^r - l) , 

Jo fM\9(u) ' 

hence by partial integration 



^ ue-™ [ X E(T y - y/p) dydx=^ (- 
Jo Jo Mw 2 V< 



\xu 2 V$(u) 
By identification of the Laplace transforms, 

Var(T*) = - f E(T y -y/n)dy. 
A* Jo 

The two inequalities in (i) now imply iii) and iv). 



- 1 



□ 



2.2. The marginal distribution under scaling. We will need the weak law of large 
numbers and an elementary renewal type theorem for T x . Such results are well-known. 
The first property below follows from the law of large numbers for X t . The second from 
the formula 

/>oo 

E{T ax )/a = / P(X at /a < x) dt -► x/n, 
Jo 

Lemma 7. As a — > oo, we have 

i) —Tax —* — i n distribution 
a /x 

ii) -E(f ax ) -> -. 
a ix 

We are now prepared to prove a limit property of the centered variable T x — xj \x under 
scaling, which is crucial for the distributional convergence in Theorem [2] 

Lemma 8. If the sequence a = a m is such that (0) holds for some c > 0, then as m — > oo, 

(17) m E{e 9(f -- ax/ » )/a ~ e 6(T "- ax ' W/o ) -> ^ f Be^^r 15 dt 

MP Jo 

and 

(18) TO ± E (e 9 ^- ax /^/ a - 1) -)■ ^- [ X e 2 e- et 'n-P dt. 

dx fi 2 /3 J 
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Proof. It is enough to prove (JT7)) since (fT5|) then follows directly from Lemma O 
Recall from d2) the relation P(X < x) = - f~ v(u, oo) dy where we use the notation 
v(y, oo) = J v(dv). For hxed cc condition on X to get 

P^ < t < f x ) = P(X a > x)P(t <f x ) + - [ P(T X _ V < t < T x )v(y, oo) dy. 

M Jo 

Multiply this identity by Be and integrate over t > to obtain 

E (e ef * - e 6T *) = P(X > x)E(e 0f * - 1) + - f" E(e ef * - e ef *-y)v(y, oo) dy. 

VJo 

Hence 

mE{e ef **/ a - e 9T -' a ) = mP(I > ax)E{e 9f "' a - 1) 

1 rx 

(19) +- / E(e eT «* /a ~e 9T ^-y> /a )amv(ay,oo)dy. 

i 1 Jo 

By©, 

— amu (ay, oo) —>■ c^y~ ~ . 

By ((U) and ©, and using the direct half of Karamata's theorem, 

mP(X a > ax) -> /3~ Vx _/5 , 

cf. Bingham ei al. (1987) Thm. 1.5.11 ii) (using in their notation f(x) — i/(x, oo), p = 
— (1 + /3), <j = 0). If we assume for the moment that the order can be interchanged in which 
we integrate over y and take the limit m,a —t oo, then applying the above asymptotic 
results as well as Lemma [7] i), 

mE{e ef -' a -e eT -' a ) 

l r P„-P(J*/H _1U l X (Jx/ii _ J(x-y)/^ r o ,-l-0 , 



Jo 

= e ««//*^ J* Be- 6t >n-Pdt, 

HP Jo 

which is the desired relation (|17p . In the remaining part of the proof we verify the va- 
lidity of this limit operation by deriving an upper bound for the integrand E(e STax l a — 
e 0T a(*-y)/ a ) ami' (ay, oo) in (fT9)l . which is <iy-integrable over (0,2:]. 
Using 

\E{e ef - - e 0f *~y)\ < \B\E[(e ef * V 1) \f x - f x _ y \] 

and Holder's inequality we have, for each integer k > 2, 

(20) \E(e ef * - e sf *-y)\ < \8\E[(e ef * V if'^ 1 )} l - 1/k E[\f x - f x _ y \ k ] l / k . 
Now, 

/"OO 

fcl 



E[\T X -T x _y\ \ - E J^ ...J 1 {T m - v <t 1 ,...,t h <f m } dt i'-- dt * 

= k\[...[ P(f x _ y <t l ,...,t k <f x )dt l ...dt 

J Ju<—<u 



>H<-<t k 

M I ... P(x - y < X tl <■■■ <X tk <x)dti... dt k . 

Iti<-<t k 
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For the event x — y < X tl < ■ ■ ■ < X tk < a; to occur it is necessary, in addition to X tl < x, 
that all increments X t . — X t ._ 1 , 2 < j < k are less than y in size. Hence the right hand 
side is at most 

k\ j ... j P(X tl <x,X tj -Xt^ <y,2<j<k)dti...dt k , 

which equals 

k\ / dt t P(X tl <x) dt 2 P{X t2 - tl <y)... dt k P{X tk ^ tk _ 1 < y) 

JO Jtx Jtk-i 



= k\E{T x )E{T y ) k 

since the increments of X(t) are independent and stationary. By (|2U|) . 
| B ( e flT«/a _ e ef a ^_ y) /a^k ^ k] ^k E [( e eT M /a v ^/(k-i)] ^^(T^/a^T^/a)* -1 . 
By Lemma [7] we may assume 

E [( e f>T a Ja v ijVCfc-i)]*- 1 ^^/,,) < 2(e te /" V l) fc (z/u), 
and thus 



A- 



for a > a and sufficiently large o , with C e ,k{x) = |0|(2 fc!) 1 / fc (e S:E ^ V 1) (x/^) 1 ^. For 
the integrand in (fl9| we have obtained 



(21) \E(e 9T °* /a - e 6T ^-y> /a )\ amv(ay, oo) 

< Ce,k(x) EiTay/a) 1 ' 1 ^ amv{ay 1 oo), < y < x, a> a Q . 

We split the further task of estimating the right hand side in the above expression in the 
two cases ay > ao and ay <a<y. 

By ([1]) and (J7J, there exist a constant C\ such that for any q < a we have ^(A) > 
Ci(AAA 9 ), A > 0. The lower bound in (JT)) ensures, moreover, that we may take q such that 
(3 < q < a. In combination with Lemma[6]i), this yields, for such q, E(T X ) < Ci{x V x q ). 
Thus, 

E{f ay /a) < C 2 (y V (y q /al' q )) <C 3 y q , 0<y<x, a> a . 

Furthermore, since the function v(x, oo) is regularly varying at infinity with index — (1+/3), 
we have for ay > ao and e > the Potter type bound 

amv(ay, oo) < C±y~ ~^ m&x(y e ,y~ e ) 

(Bingham et al. (1987), Ch. 1.5). Thus, for some constant C, 

E(f ay / a) 1 ' 1 /" amv(ay, oo) < C y^ 1 ' 1 ^ y~ l ^-\ 

Since q > f3 we may take k so large that q(l — 1/fc) > j3 and then e so small that 
e < q(l — 1/fc) — /3 to obtain a dominating function for the integrand in (|21[) which is 
integrable in y over [0,x]. 

For the remaining case ay < ao, Lemma |6]i) implies 

S(f 02/ /a) 1 - 1/fe am^(at/,oo) < {e 2 /{e - l)) 1 ' 1 ^ ma 1/k y{ay, oo) x ' k . 

Using a property of slowly varying functions (Bingham et al. (1987), Prop 1.3.6), for any 
e > 0, L(a)a e -> oo as a -> oo. Hence we may assume a~ e < L{a). Also, 



v{ay, oo) < — / uv(du) < — . 
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Thus, using ([B]), 

ma l ' k v{ay,oo) x ' k < I^) *«H-i/* ( M / a j,)V«= < 2cV +1/ *a e+ V 1/fc - 
Now apply a < a^/y to obtain from (|2ip a constant C for which 



|£:( e eT aa: /a _ e eTa(»-„)/«)| aml/ ( ayj0O ) < C a% 



+0 y -P-l/k-e^ 



This is again integrable if we make the same choise of k and e as above. This concludes 
the proof that the limit in (JTHJ) can be carried out under the integral sign and hence the 
proof of the lemma. 

□ 

Lemma 9. For a — a m such that (0) holds for some c > 0, 

mE ( e 6{T alc -ax/n)/a _ ]_) _> __ / / p^Bs/ » s -fl dgdy ^ m _^ OQ _ 



P^ 2 JO JO 

Proof. By Lemma [SJ mE{e 9 ^ Ta:c ~ ax ^^ a — 1) is nonnegative and increasing in x. The 
limit function on the right hand side is also nonnegative and increasing. Hence the lemma 
follows from weak convergence of measures if we can prove 

(22) / e - ux ± m E(e 0{T -- ax/ ^ /a -l)dx-> / e- ux ^— \ (Fe^'^a^ ds I dx. 
Jo d x Jo VP < P Jo ) 

To find the Laplace transform on the right hand side note that 

« 2 r(i-/3) e 2 f°° _ ux Bx a 

' e ux e ax x p dx, 9 < u. 



Piu-e) 1 -? p J 

Multiplication of the transform by 1/u corresponds to integration of e 6x x~^ . Hence 

f°° -ux (}_ r q2 0s -0 A d = T { 1 ~ P)° Q <u 

and hence (|2"2"|) is equivalent to 

(23) ^ rt£(e 9 ^-W/" - 1) <fe -► T ( {1 ~ P K*f> 2 > ° > -/-• 

Jo (3u(u + 9/fi) 1 Pfi- 

To help analyze the Laplace transform in (|23|) we introduce the additional notation 
I(u) = (M - $(«) = / (e' ux -1 + ux) v{dx) > 0. 



Jo 

,2 r°° D -ni7 



Writing I(u) — u 2 L e ux U(x) dx with U(x) — J v(y, oo) dy, it follows from Kara- 
mata's Tauberian Theorem (Thm. 1.7.6 in Bingham et al. (1987)) that 

T , , s T(l - P)L(a/u)u 1+l3 

(24) al(u/a) 1 'J , a -> oo 

pa p 

Relation (fblj) of Lemma [4] now shows 

ue- ux mE{e e{T ^- ax /^/ a - 1) dx 
m9 2 al({u + 9 //x) /a) 



(fiu + 9) 2 fm - al((u + 9/fi)/a) 

mL(a) (1 _^) T{l-P)9 2 

1 — (u + On) — 5-3 

a p pn u 

r(i_/3) c /3fl2 



puiu + o/n) 1 -^' 



X _P 2 , > -fiU, 
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which proves (|23[) and hence the lemma. 

We are now able to conclude convergence of the marginal distributions. 

Lemma 10. Under the assumptions of Theorem, 2, for any x > 

m 
— ^2( T a m x - -OmX) -> --cYp(x/c) 

Proof. Writing 



□ 



;=i 



A im \0;x) = m E(e 0{T ^- ax/ » )/a - 1), 



Lemma [9] shows that 

i m 1 

log £ exp {e— J2( T a k J x ~ -a m x)\ dx 



fc=l 



log(l + — A( m H8;x)) m ^ 

V 777, / 



,/3 i-x c-y 



■)2„-e s /n Q -p 



m > Pii 2 Jo Jo 

This proves the lemma since the limit process Yp has the property 

1 f x f y 

9 2 e 9s s- p dsdy 



s p dsdy. 



\ogE(e eY ^ x) ) 
and so, as noticed in Lemma [21 

\ogE(e- ecY ^ x / c ^n = ^ [ X I" 8 2 e~ es ^s-P dsdy. 
PH Jo Jo 



P Jo Jo 

r f3 fx ry 



D 



3. Multivariate distributions 



The proofs of convergence of the finite-dimensional distributions are based on the fol- 
lowing recursive equations for moment generating functions. 

Lemma 11. Fix n > 2 and a sequence of time points < X\ < • • • < X n . The moment 
generating function of the finite- dimensional distributions of the stationary inverse Levy 
subordinator process {T x } satisfies the recurrence relation 



n n 

E exp { Y, OiT Xi } - E exp { ^ Q{T Xi } 



(25) 



E 



f Xl r r " ~ ~n r r " ' i 
j^-j- / E [ exp | J2 6iT Xi - x }J d x E [ exp |T B J^ 9 i\ 



i=2 i=l 

where T x is the corresponding pure inverse Levy process. Moreover, 



11 71 

E exp { J2 e iTxi } = E ex P { I] fcr B1 } 



4=2 



(26) 



E 



^— / E [ exp { J2 9if Xi _ x }] 4£ [ exp {t x ^ 0*} 



Proof. We have 

n n 

E exp { J2 QiTxi } - E exp { ^ fcT., } 

i=l i=2 

Since 



exp | / y i/a- , 

i=2 



^ ^ -1 



e x i — 1 = 



1 



{n<T xl }Wie 



du= I l{x u < Xl }0ie 9lU du, 
Jo 
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i:s 



it follows that 
E 



j=2 

r r°° ( n i 

n n 

l { x u <, l} exp { J2 e *( T ^ -Tx u )} <?i exp {(£ ^)T^} 



d» 



i=2 



Here, Tx u = u. For any u > and i > 2, on the set {X u < x±} we have 

{T Xi - T Xu < *} = {T Ki <u + t} = {X u+t > Xi }. 

Since {X t } has independent increments the rightmost event has the same probability as 

{X u +Xt>Xi} = {T Xi _ Xu <t}, 

where X u < x\ is assumed independent of X t . Thus, on {X u < x\\ the increment 
T Xi — Tx u has the same distribution as T Xi ^x u - It follows that 

n n 

E exp { J^ OiT Xi } - E exp { ^ d(T Xi } 

*=1 i=2 

/* oo _ ^ ^ 

= M?[ / l { x u <. l} £[exp { J2 OiT xt - X J\\xJ exp {(^ ft)«} 



4=2 



(/(/ 



(=1 



q x e 



/ E [ exp { £ ftf 01 _ a }] exp { (^ Z )T X } dT x 



i=2 i=l 

where the integration after variable substitution x = X u is with respect to the increasing 
function of bounded variation {T x , x > 0}. (Intuitively, the time-change X u picks out the 
rightmost point of each flat piece of T x .) Moreover, if we change to the measure 

n n n 

d x ( exp \t x J2 0i}) = ( 12 d ) ex P { T * 12 9 *} dT * 

2 — 1 l— 1 i— 1 

we obtain 

n 



E exp { Y, 0iT Xi } - E exp { J] 0*7^ } 



;=i 



£ 



i=2 



En 
i=l 

"ECU £ Mg«"} 

which is (|23|) . Start with T rather than T to get 
For n > 1 and 1 < fc < n, put #£,„ = (0^, 
= xq < x\ < • • • < x n and let 



I E [ exp { J2 9iT* t - x }] d x ( exp \t x ^ Z }) 



n 

{r*5>} 



daB 



exp 



□ 



9 n ) and x fc ,„ = (a;*., . . . ,i„), where 



(27) 



$, 



-fe+l(."fe,n! x k, 



it 

1 )= J Bexp{^^(T Xi -x,/M)} 



denote the multivariate moment generating functions for the centered process {T x — 
x/ (i} x >q. Here, the subindex n — k + 1 is the number of elements of the argument vectors 
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@k,n> %k,n- Similarly, let Q n -k+i(0k,n', Xk.n), 1 < k < n, denote the corresponding func- 
tions for the pure process {T x — xj fi} x >o- The subtraction Xk, n — u = (xk — u,...,X n — u) 
is interpreted component-wise in the next statement and in the sequel. 



Lemma 12. The moment generating functions defined in j[27\} satisfy the integral equation 

*»(01,»5 Sl,») = *n-l(ft»,n; ^K^ 1 ^ 



2_/i=l &i Jo 



-Ql(xi—x)/n 



<&n-l(02,n; X2,n - X 



71 

)*i(£«<5 

1 / e- fll ( Xl - x )/"*„_i(ft ! , n ;S2,n-a:)*i(y)ft;a!)da: 

4 •'o \ =1 7 



M Jo 
Proof. By Lemma fTTl 

$»(0i,»; xi, n ) = $„_i(0 2 ,»; ^e-* 1 * 1 ^ + 



sr=i A* 



< / e- e ^- x ^^ n ^ 1 (e 2 . n ;x 2 . n -x)exp{--Y / S l }d x E[exp{T x Y,^} 
which, by observing 

n n 

exp {-- £ 4 } d x ^[exp {r x £ft}] 

n -. n n 

= d x #[exp |(T X - x/m) £ 9i}] + ~ J2 ^[ CX P {( T - " x/fi) J2 9 *} 

i=l " i=l " i=l 

n 1 n n 



d.v 



i=\ i=\ 



n 



may be rewritten in the form stated in the lemma. 

According to (JT2J) we must find the limits of the scaled function 

n 

m($ n (§i t n/b; axi^n) - 1) = mE exp | £ 9i(T aXi - ax l /^)/b\ - 1 

as m, a and b tend to infinity, when a and & satisfy either ([5]) together with ((5J) or condition 
([5]). The first case is FBM scaling leading to fractional Brownian motion in the limit, as 
in Theorem 1, and the second case (with a — - b) is the intermediate scaling studied in 
Theorem 2. 

For n > 1, m > 1 and a, b > we introduce 



as well as 

and 

(28) 



Ak m) (0i, n ; xi, n ) = m($ n (9i tn /b; axi, n ) - 1), 



A< m) («i,„;ii,„) = ^-{®n{0i,n/b;ax x , n )-l) 
b 



Si m) (£i,n; *i,n) = A( [ m )(0 1 , n ; xr,„) - A^&y, z x , re ). 



Our strategy for finding the corresponding limit functions is to derive for fixed m sequences 
of integral equations, which are recursive in n. As already pointed out we give the detailed 
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proof only for Theorem 2. To simplify notation during this analysis we will also use a 
special notation for the sums 

n 

Vn =X^ 4 ' U ~ 1 ' 
i=l 

which appear frequently as evident in Lemma [T2l 

3.1. Multivariate distribution under the intermediate scaling. We study the as- 
ymptotic limits of An and A„ as m — > 00 under assumption (|6]). For simplicity the 
constant in ([6]) is set to c = 1. The general case c^ 1 then follows from Lemma [2] We 
begin with a system of equations for the functions S„ ' defined in (f2"5|) , which will be used 
to determine corresponding limit functions as m — > oo. 

Lemma 13. We have 

HOi/Vn - 1) r e-O^-^E^ivn, dx) + -R^ m \ 
Jo m 

where 

R {m) = - [ X \- e ^- x)/ »At\(kn,X2,n-x)Zt\r l n, dx) 
Vn Jo 

Proof. After inserting the scaling parameters m and a into the equation obtained in 
Lemma [T^] and sorting the terms appropriately, it is seen that the scaled functions A„ 
and A„ satisfy 

\( m )t5 ■ * ^ — a(" 1 ) la ■ * \ c S lXl /u , r(m) T {m) T {m) 1 p (m) 

where 

% Jo 

4 m) = - r e- e ^-^^At\(hn;x 2 , n -x)dx 
M Jo 

4 m) = ^ P e' 6 ^-^^K { r\vn\x)dx 
M Jo 

4 m) = - f 1 e^'" 1 -)/" A^fon; x 2 ,„ - or) A< m) fa„; dx) 

Vn JO 

9 1 r^fxi-^MlWM __.= _ _,UW/ 



By Lemma [5j 



e-" 1 ^-^/" An(fl 2 ,„; x 2 ,„ - s) Ai"'^,,; x) dx 
M Jo 



4a("°(0; x) = ± m E(e e V«- m 'M a _ !) 

ax dx 



(29) = ^ TO £( e e ^- a ^)/ Q - e e(r„-»/*«)/a) = ^~( m )(0 ; x ). 

Thus, 

4™) = ^ r e-flxCxx-.:)//. A^\(fl a , n ; x 2 ,„ - x) A[ m \ Vn ; x) dx. 
f 1 Jo 
By partial integration, 

/M = A^fo; an) - f * c -*0«- ^A^fo; dx). 
Jo 



16 INGEMAR KAJ AND ANDERS MARTIN-LOF 

Hence 

Ai m) (0"i,„; x x , n ) = ?Sr\{hn, X2,n)e- e ^ /f * + A { " l) (Vn; si) 

(30) +4 m) +4 m) + -4 m) , 



77) 



where now 



Similarly, 



(m) 



(^--l) f 1 e- 6 '(^)/"Af)( I)n ;&). 

V ??rc I Jo 



(3i) +4 m) + 7i m) + -4 m) , 



with 



and 



r (m) 



(^-i) r e -M«i-.)/4w^ ;(b ) 



XI 



4 m) = — / e - 9l ^-»/"Ai , ! ) 1 (e 2 , n ;^ !n -.)AW(7 )n ;&) 

»7« Jo 

+ * r\- e ^^y^At\(e2,n 1 X 2 . n ~x)A ( r\vn; x) dx. 

M Jo 

By subtracting (|30|) from ([31]) and using i?( m ' = R^ 1 — i?-f m we obtain the desired 
equation for S„ . 



□ 



The next result generalizes Lemma [5] to the multivariate distributions. 
Lemma 14. For any < s < %\ and n > 1, 

--^-Al m )(0 1 ,„; x ljn - a) = ^4 m) («i.n; Si.„ - s). 
as /u 

Proof. For n = 1 this is P5]). For n > 2 and < s < xi, using (j50|) . 
-f A™(0 lin ; x 1;n a) = ±±R[ m \ Xltn a) + -^-A{ m) (vn, *x - *) 

LAO //i LAO '/?7, U6 

+ (fi _i)^l [ X1 - e -B l{xi -s- x )/M A (-) (77n . da.) 

V ??n ' JU Jo 

+e-^^-)/"^A^\(fe,„;« a , n - a) - ^ e -*C»i- )/a- S^^;^ - a) 
as fi 



and hence by ((29 



~A( m )(^, n ; S 1;n - a) = --^-Rt\x 1>n - a) + ^S^fa, a* - 
as m as [i 

-e- fl1 ^-)/" (^A^l^a.njxa.n - «) - ^H^x^;^ - s) 
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Here, 



x± — s 



vi e °^z ( r\ri n ;x)dx 

M Jo 



e ^(-x-)/M S (-) ( ^. Xl - S )- " e^H^fo; dx) 



X\—S 



(I 



S( . 



ds to ds |U 

-fo»-0i) f "e-°^-°-'V'*Ei m \ Vn ; dx) 



V Jo 



,-e^-s)/, | * A£\(0 2 , n ;5 2 , B - «) - ^i!'!^,;^,,, - a) 
\ds /i 



By replacing the integral term using Lemma [T3] this implies 

. 1 n 

-±Ai m )(^, n ; £ lifl - a) = -5]^4 m) (^i, n ; «i,n " « 

2 — 1 



l(Kl " S)Ai (l A i- ) i(^,„;S2,n - a) + ^^si- } ifen;x 2) „ - a)\ 

1 1 ™ 

i£ m \xi, n -s) — J2 eiR{m) (^> n - s )- 



3--R1 (*i,n -a) 

m as to /i 

Since 

M m) (*i,n " a) = - f 1 e- 9 ^-^At\(e2,n,x 2 , n - x)A[ m) ( Vn ; x-s)dx 
A* J s 

we have the identity 

-— R\ m '(x 1}n - S) 

= °1 l^ e- e ^-^^At\(hn;X2,n ~ x)d x A ( r\vn, X s) 

M Js 



1 " 

= -J2eiR( m \xi, n - S ). 

Thus, 



**<=i 



-A^\h. n \ X!, n - a) + ^ 4 m) (^i.«; &!,„ - 5) 
as /i 



e -6i(xi-s)/n 



( ^-At\(02, n ; x 2 . n - s ) + ij2^ zt\(kn;x2. 



The statement of the lemma now follows by induction. □ 

Lemma 15. For each n > 1, the limit functions 

S„(0i,„; x 1>n ) = lim S^ m) (0i jn ; zi i7l ) 
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exist and are given by 

n j—X 
H n (01,„; Xl tn ) = ^ ex P { - y^(fli H 1" 0n)(£j - ^i-l)/^} 

(32) x - f ° {9 j + • • • + 0„)e- ( ^+- +e " )(u -^- l)/A '/3- 1 w- /3 du. 

A* Jxj-i 

For n > 2 i/iej/ sofce i/ie recursive system 

S„(0i jri ; xi jn ) -Hi(r?„; xi) = e~ eiXl/ WH„_i(6>2,n;S2,70 -Hi(?7 n - 9\\ xi)V 
Proof. As m — ^ oo, by Lemmas [5] and [51 

(33) H< m) (#; x) - !iJL m E(e e ( T »- a *'M a - 1) -> -J- f fcr^/"^ du = H x (0 

6/ da; P)U J 

Using (|2"5)) . for arbitrary a, 



o 

= Z[ m \a; x x ) - ^ r e- e ^-^^E[ m \a; x) dx 
M Jo 

a Jo 
= S< m) (a; xO - ^A[ m \a; xi) + ^- H e -*(»i-*)//' A< m) (a; x) dx. 

Here, (|33[) shows that H^ converges to Hi and Lemma [5] shows that A^ m ' converges to 

a limit function Ai which satisfies 4-Ai(a,x) — aHi(a,x)//i. Since < A" 1 (a; x) < 

Aj" 1 (a; Xi) for < x < xi by Lemma^ we can find a dominating function for Aj m (a; x) 
on [0, xi] and conclude that the last integral term also converges. By reverting the partial 
integrations this yields with a = r\ n 

(34) r e-'^-'VuE^fa, dx) -> H e^^^Hi^; dx). 

Jo Jo 

Let us now consider the remainder terms i?( m ) in Lemma Q2J Because of (j2"5)l we may 
rewrite using 

A< m) (a;dx) - A ( ™\a;dx)+z[ m \a;dx) 

= A { ™\a; dx) +e- ax ^d x (e ax ^ l E^ ) (a; x)) - -S^a; x) dx 

M 

= e- ax/ ^d x {e ax ^E < ( n \a;x)), 
and obtain, again with a = J2i=i @i-> 

R {m) = 9 j-e-^/» p e-^-^^At\(02y,x 2 , n -x)d x (e a ^S^(a;x)). 

The above integration is carried out with respect to the function 

F^ m \x) = e"/"a[ m) (aj x) = mE(e a¥a * /a - e aT - /a ). 

It was observed in the proof of Lemma \5\ that T can be viewed as a shift of T with the 
first point of increase in Xq. In particular, T x — T x equals Tx almost surely on the set 
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{X < x}. Hence the increment 

£( e afWh _ e <*T*+K) _ E(e afx - e aT *) 

= E{{e aT "+ h - e aT *)(e aTx o -l),X <x) 

+E(e aT *+ h (e aTx « - 1) - (e aT * - 1), x < X < x + h) 
+E(e aT *+ h - e aT *, x + h<X Q ) 

is positive for a > and negative for a < 0. Thus, F^ m > is a monotone measure with limit 
e Qa: ' M Si(a; x). Using the variation measure i 7 '^ r "'^ | we obtain a constant C n {x\) uniform 
in m such that 



\R {m ^\ < 



sup 

0<a;<a:i 



lX1+{a - ei)x) ^A^\{6 2tn ;x 2 , n -x) \F^\( Xl ) 



;m 



< C„(xi) SUp \M n _ 1 {9 2 ,n]X 2 ,n- x)\ 
0<x<a:i 

The same arguments apply to the remainder terms R™ in (|30[) and R 2 in ([31]) . Hence, 
taking C n {xx) sufficiently large, 

^ m) <C„(xi) sup \A { ™\(e 2 , n ;x 2 , n -x)\, i = l,2. 

0<x<xi 

We are now prepared to carry out an induction on n in Lemma fTS] and equation pip . 



77M 



'Am) 



Assume that 'E n _ 1 {9 2 , n ]X 2: n) converges to S„_i(# 2 ,n; x 2>n ) and X n _ l {9 2 ^x 2 ^ n ) is such 
that 

SUp |A^\(^ 2 ,n;»2,n ~ »)! < Cn-l(&l,n)- 
0<ie<:ei 

Then, applying (l3"3")l . (jM|) and the induction hypothesis to Lemma 02] it follows that all 
limit functions 2 n exist and satisfy 

3 n (0i,n; Si,„) -3i(t7„; Xi) = e^^^Hn-i^.nS^.n) 

+ (^-l) r e -^(— )/" Sl (^;(fa) 



r(m) 



which is the desired relation. Moreover, observing that the convergence of -Q is a 
byproduct of the proof of ([34]) , it follows from (|3Tj) that 



sup \A^(e 1>n ;x 1<n -u)\ <C n (x , n ), 0< xq <x x < ■■■ < 



To verify the explicit form (1321) of the solution, assume that the claim is correct for 
index n — 1. Then 



/"Jo 

n-l j-1 

^expj - (6 2 -\ h 9 n )x 2 /n - ^(6» i+ i H \- 9 n )(x l+1 - x. t )/[i\ 



i=2 

x — 



i=1 



J 3 + 



! + ••• + 6! rl )e-^'+ 1 +- +e ">( u - :E ^ t /r 1 u-' 3 dw. 
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This implies 



4=2 

i r x 2 



i=2 

9 2 + • ■ • + ^)e" (e2+ - +e " WM /3~ 1 w-' 3 du 

i-i 
+ £ exp | - (0 2 + • • • + 0„> 2 //" - 5^(0i + ' " " + *»)(*' " aJi-0/A*} 

4=3 

+ ■ • • + n )c-( fl i+-+ fl »)( tt -»i-O//*^-l„-/» ^ 



j=3 4=3 

1 /"^ 



Hence 



9ia;i/ ' 1 (s n _ 1 (0 2 ,„;5 2 , n )-H 1 (^;^; : 

4=2 

= e -(«i+-+fl„)»i/Mi r 2 ( fe + . . . + n ) e -(02+-+0„x«-*i)//^-i u -/3 dw 

n J — 1 

+ ^exp| -^(0, H \-0 n )(xi -Xi-i)/n\ 



x- / (^■ + --- + 6 , „)e- ( ^ + - +e " )(u -^- l)/ ' i /3- 1 u-^du 



j=3 4=1 

1 /^j 
n j—\ 

^cxpj -^2(9i H h6» n )(xi -Xj_i)///j 

3=2 4=1 

1 /"^ 



x — 



A* Jxj 



+ ■■■ + 0„) e -(^+-+ e ")("-^-i)/"/3- 1 w -/ 3 d u . 



n 
The remaining proofs of the convergence of multivariate distributions in Theorem 2 are 
organized in three consecutive lemmas, leading up to the identification of the cumulant 
generating function dill) in Theorem 2. 

Lemma 16. The limit functions A„ = lim TO _ s . 00 A„ , n > 1, exist and we have 

nl rx i-u 

(35) K 1 {9-x) = -—r / e-^'^v-P dvdu 

P> Jo Jo 

and for n > 2, recursively 

A„(0~i,„; 5i,„) = A n _i(0 2 ,n;S 2 ,„ -xi) + Ai(?7„; xi) 

+ (— " P e-^^-^^Ar^; ete) 

V ?7n ' Jo 

-(— 6± -h- + l) / Xl e- 91 ^ 1 -^^ dxA n _i(0 2) „; a: 2 , n - x). 
\r)n-(>i 'Jo 
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Proof. For n = 1 this follows from Lemma H] and for n > 2 from (|30p and a further 
partial integration of the term 1% , which gives 



Ai m) (#i,„; xi, n ) = ASd.; 5 2 ,„ - a*) + A^V; x x ) 



TJn ' Jo 

:r i 

e 

ii 

Xl 



^"'^iAfife™;^,-!) 



H J to 

= A^™'^^,™; x 2 , n - Xi) + A\ m \r] n ; Xi) 

+ (_! - A r e - e ^-*V»A[ m \r) n ; dx) 
V??« / Jo 

-(— ^V + l) r e- eiixi -* W d x At\(02,n;Si2,n - x) + -R^\ 
\rin-9i 'Jo m 

where we apply Lemma Q3] for the last equality. The arguments which justify that we are 
allowed to exchange the order of integration and taking limits in to, as well as controlling 
the remainder terms, are parallel to those in the proof of Lemma 1151 again based on 
Lemma [Til □ 

In view of (|12p we conclude from Lemma [16] the convergence of the finite-dimensional 
distributions in Theorem [2] 

Lemma 17. The finite- dimensional distributions of the sequence of random processes 
studied in Theorem^ (with c = I) converge to those of a limit process Yp, such that the 
collection of logarithmic moment generating functions 

n 

A„(0i, n ; Si,„) =log£ , exp|^6» i y" / j(x. i )|, n>l 
is the unique solution to the closed system of linear integral equations 

n 
A„(01, n ; Xl, n ) = A„_i(6> 2 ,n;£2,n ~ *l) + Al ( ^ dj-, XiJ 

i=l 

+ ( *L_-l) re-^-^AJ^dx) 

K l^i=l °t y Jo K i=1 ' 

-(— *l_ + l) r e- e ^-^^d x A n ^(e 2tn] X2 t n-x), n>2, 

V 2^=2 ti i ' JO 

with Ai as in i35\) . 

Lemma 18. The cumulant function for the increments ofYp, 

n 

r„(0i,„; xi, n ) = log^exp|^6' 4 (y / 3(a; 4 ) -Yp(xi-i))\ 

i=i 

has the explicit form given in ill]) . 
Proof. We have 

rVj(#l,n; Xl,n) — A n {{0i — #2, ■ ■ ■ , &n-l ~ Q n ,Qn)i %l,n) 
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so by Lemma [T71 

r„(£i,„; xi,„) = r„_i(0 2 ,n; S 2 ,„ - xi) + Ai(6»i; xi) 

-Pi-feX*!-*)//*^^; (fa) 
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7 i Jo 

e-^-^^-'^dTn-i^njsa,,, -a:), n > 2. 

7 2 Jo 

It may now be checked that the functions in (fTTj) solve the above system of equations. For 
details, see Gaigalas, Kaj [5], Section 6.3. □ 

4. Remaining proofs 

4.1. Limiting distribution under FBM scaling. In this section we discuss briefly the 
convergence of the finite-dimensional distributions in Theorem 1. Recall that for standard 
fractional Brownian motion B^ , 

n _. n n 

log £/ exp | ^ ^Bh^) | = -cr 2 , ^353^ Cov{B H {x i ),B H (x J )) 

2—1 2= 1 J — 1 

with 

Cov{B H {x),B H {y))^^{x 2 -f i +y 2 -f i -{x-y) 2 -' 3 ). 

In the scaling regime defined by ([5]) and © we have 



a / a^fi am a^m/J, 

By analyzing in this case the recursive equations for An (@i,n] %\,r) it follows that 

r2 n 

".J 



(36) A n (0 li „;S 1 .„)= lim A% n) (6 h n;x 1 , r ) = ^y2 6 i xl i 



■P 



nh ' %i ' ""-fti-tir 



i=\ 



Moreover, the limit functions A n (9i. n ;xi. r ) = lim m _>. 00 A„ (^i,n5*i,r) satisfy, in analogy 
to the result of Lemma [T71 



Ai(0; s) = lim A^ m) (#; x) = \alp,- 2 B 1 x 1 - 



and for n > 2, 



A„(6>i,„; xi )tl ) = A„_i(6> 2 ,„; x 2 ,„) + -^w n Ai( y^fl»; x\ 






9 f Xl ~ 

— A n _i(6>2,„;£2,« -x)dx 
M Jo 

Thus, using (j3"6"j) . 



n- 2 n i 

°s v-^ „ „ i r 






Hence 

2 rc n 



a . * 1 
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4.2. Proof of tightness in C. To complete the proofs of our results we establish tightness 
of the sequences 

1 m 

y<™> (a . ) = J_£(rW,-^) 

studied in Theorems 1 and 2, by applying a standard moment criterion. Since Y^ m \x) 
has stationary increments, to prove that {F^™)} is tight in C it is enough to find 7 > 1, 
an integer too and a constant K such that for fixed T, 

m 



(37) Var(F( m )(x)) = ^-Var(T amX ) < Kx 1 



for < x < T and to > to (Billingsley (1968), Thm. 12.3). 

By Lemma [6] iii) , the variance of T x is a non-decreasing function in x. Hence we may 
apply Karamata's Tauberian theorem (Bingham et al. [4] Theorem 1.7.1) to show that 
Vax(T x ) is regularly varying in infinity with index 2-/3. Indeed, recalling the previously 
used notation I{u) = /au — $(u), the asymptotic property (j24l) implies 



^■)-iE©"~"^^ -0, 



, (/iu) 2 ^ V ^u / /3n 2 u 2 P 

hence 

(38) Var(T x )~ r(3-/3)/3 M 3 " M 3 ( } ' 

The next step is to apply the Potter bounds for regularly varying functions (Bingham et 
a/., Ch 1.5) to obtain for any e > an ao, such that 

Var(T '" ) < (1 + e) maxOr 2 -^, x 2 -^), a > a , ax > a . 
Var(T a ) 

Hence for to > too so large that a > ao, ax > ao, 

mV&r(T ax )/b 2 < (I + e)mVar(T a )b- 2 max(x 2 - 0+e ,x 2 - p - e ) 

But in either case of the FBM scaling ([5]), © or the intermediate scaling ([5]), the asymp- 
totic relation (|3"5|) yields 

mV&l(T a )/b — > (7a/ '/i , to — > 00, 

and so, eventually choosing a larger mi > too, 

TOVar(T aa; )/6 2 < (1 + e)(cr^ + e) max(x 2_ ' 3+e , a; 2- ^ -6 ), to > mi. 



With e < 1 - p this yields ([37]) for as > a . 

It remains to prove (|37p for a > a and ax < ag . By Lemma |5] iv) , 



^ />ax 

mV^(T ax )/b 2 < -g / $(1/2,) 



n , 2e to ax 



ix b 2 ®(l/ax) 

By (Jll) we can find a constant Ci and q > j3, such that <j>(\/ax) < Ci(ax) q . As in the proof 
of LemmaEJ we may take a _c < L(a). In Theorem 1, ma 2 ~^L(a)/iib 2 = 1. In Theorem 
2, ma 2 -^L(a)/fxb 2 -4- c' 3 . Thus, 

Since (3 < q < a <1, we may take e < 1 — q to obtain (|3T|) for 7 = 1 + g — /3 > 1. 

□ 
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